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^ We have measured the persistent current in individual normal metal rings over a wide range of 

CN magnetic fields. From this data, we extract the first six cumulants of the single-ring persistent 

I — ^ current distribution. Our results are consistent with the theoretical prediction that this distribution 

^ should be nearly Gaussian (i.e., that these cumulants should be nearly zero) for diffusive metallic 

c/3 rings. This measurement highlights the particular sensitivity of persistent current to the 

^ mesoscopic fluctuations within a single coherent volume. 
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One of the fundamental archetypes in mesoscopic physics is a system whose volume is 
sufficiently small that electrons remain quantum coherent within it, yet sufficiently large and 
complex that its energy spectrum cannot, in practice, be calculated exactly. Such a coherent 
volume can be realized, for example, in a microfabricated metal device cooled to sufficiently low 
temperatures. In such a device the electronic spectrum will depend upon the angstrom-scale 
disorder within the metal, which is beyond the control of most fabrication techniques. As a 
result, the device's spectrum (and all related physical quantities) will be drawn randomly from an 
ensemble representing the possible realizations of the disorder within nominally identical (i.e., 
lithographically identical) devices. 

The sample-to-sample fluctuations that result from this effective randomness are 
characterized by a distribution function Px. Here x represents any physical quantity that depends 
upon the disorder, for example the dimensionless conductance g or the persistent current /. These 
distributions play an important role in our understanding of how electrons flow through 
disordered materials, for example in Anderson localization and the scaling theory of 
conductance. ^'^'^ For metallic samples (i.e., with mean conductance (g) » 1) calculations 

predict that Pg and Pi approach a Gaussian distribution as (g) ^oo ."^'^'^'^'^'^ Deviations from 

Gaussianity at finite (but large) (g) reflect the approach of Anderson localization and the 

"breakdown" of single-parameter scaling that is due to a finite-sized system's vestigial 
sensitivity to the particular details of its microscopic disorder.^'"* 

Several measurements of Pg, the distribution of mesoscopic conductance fluctuations, have 
been made in systems with (g) < 1, including ballistic semiconductor quantum dots^*''"'^^ and 

semiconductor wires near the localization threshold. In this regime, experiments have found 
agreement with theory. However in metallic samples, experiments to date have largely been 

confined to measurements of ((g^ )) , the second cumulant of Pg. These measurements of Ug]) 



have found excellent agreement with theory in a broad range of circumstances. > ' ' ' 
However, little is known experimentally about mesoscopic fluctuations (of g or any other 
quantity) in metals beyond the second cumulant. 

Measuring the full distribution of mesoscopic fluctuations in a metal device is challenging. 
In part this is because most experiments detect g, and so attach leads to the device. These leads 
are much larger than the electrons' phase coherence length L^, and so contain a large number of 





coherent volumes that contribute in some degree^'^ '^'^*^ to the measured g. (The contribution from 
the leads is less important when g ^ 1, as in Refs.[10,ll, 12,13].) Since the fluctuations of each 
of these coherent volumes are assumed to be independent, higher cumulants of Pg will tend to be 
suppressed in such a measurement, with the result that the observed fluctuations will appear 
more Gaussian than the actual Pg. The impact of the leads can be reduced by measuring Pg in 
long wires (i.e., much longer than L^,) but this ensures that the wire itself contains many coherent 
volumes, with the result that the observed fluctuations will again appear more Gaussian than 
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However it is possible to measure the mesoscopic fluctuations of a single coherent volume 
by detecting the persistent current / in an isolated metal ring. This has been challenging in the 
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past owing to the small signals involved, but it was recently shown that micromechanical 
torsional magnetometers can measure persistent current with very high sensitivity and low back- 
action.^^'^^ This technique has been applied primarily to arrays of metal rings, with the result that 

the first two cumulants ^(/)^ and ^(^^)^ of Pi were measured with high precision over a wide 
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range of parameters. ' However the sensitivity achieved in Ref. [23] (as well as in other studies 
of individual metal rings^^'^^) did not allow for measurements of individual rings with adequate 
signal-to-noise ratio (SNR) to resolve the higher cumulants of P/. 

Here we describe measurements of the persistent current in a large number of individual 
rings. From these measurements we extract the first slk cumulants of Pi, as well as other higher- 
order statistical properties of the persistent current. This is achieved by improving the SNR of the 
technique described in Refs. [22,23], and by combining data from more than 400 effectively 
independent measurements. We find that our results agree with theoretical predictions to within 
the sensitivity of the measurements. Specifically, we find that that the first six cumulants of Pi 
are consistent with a Gaussian distribution. The small deviations from Gaussianity that are 
predicted by theory are too small to be detected in our experiment. 

A tj^ical device is shown in Fig. 1. It consists of a single-crystal Si cantilever supporting a 
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single Al ring. The fabrication of these devices has been described previously. ' The rings 
were fabricated via standard electron-beam lithography and were evaporated from a 99.999% 
purity Al source onto a Si substrate with a native oxide. In addition to the ring samples, Al wires 
and contact pads were co-deposited onto the same wafer to allow for transport characterization of 
the metal. Details of these transport measurements are given in the Supplemental Material (SM). 



These measurements show that L(p > Ijrr (r is the rings' mean radius) for the temperatures at 
which the persistent current is measured. They also provide the electrons' diffusion coefficient D 
= 0.020 + 0.0015 m^/s. 

The procedure for measuring the persistent current (PC) has also been described 
previously. The cantilever's displacement is monitored by a laser interferometer. The signal 
from the interferometer is used to drive the cantilever in a phase-locked loop, allowing the 
cantilever's resonance frequency co^ to be monitored. In the presence of an applied magnetic 
field B, the persistent current / circulating in the ring produces a torque on the cantilever. This 
torque changes cOm, and / is inferred from this change. Details are given in the SM and Ref. [23]. 

Measurements of the PC were made at temperatures 320 mK <T< 365 mK, and magnetic 
fields 4 T < 5 < 9 T (applied at an angle a = 45° relative to the rings' plane). This is well above 
the critical field of Al, ensuring that the rings are in their normal state. The large B is required to 
produce a detectable torque. It also simplifies the data analysis, as large B within the metal of the 
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ring strongly suppresses the effect of electron-electron interactions on the PC, allowing us to 
compare our results to independent-electron theory (though we note that interactions are not 
predicted to make Pj non-Gaussian for large (g) ^'^'^). 

Measurements were made on eight different rings, with each ring on a separate cantilever. 
The full data sets from each ring, as well as the rings' dimensions, are shown in the SM. A 
typical measurement of I{B) for one of these devices (device #6) is shown in Fig. 2. The rapid 
oscillations in Fig. 2(a) & 2(b) are due to the Aharonov-Bohm (AB) effect: as B is varied, the 
magnetic flux O = BA^sina through the ring varies, causing 1(B) to oscillate with period 5per = 
Oo/(A^)Sina). Here A,d = is the typical area enclosed by the electrons in the ring, and Oq = h/e. 
No higher harmonics of the AB oscillations were observed above the noise floor of the 
measurement. 

The AB oscillations' amplitude (and upon closer inspection, their phase) varies on a field 
scale larger than 5per. These aperiodic modulations result from the fact that sweeping B also 
varies the amount of magnetic flux in the metal of the ring <I»in ~ iTurwB, where w is the ring's 
width. Because the ring represents a single coherent volume, its spectrum is expected to be 
randomized each time Om changes by ~ Oq. The ergodic hypothesis identifies this 
randomization with a new realization of the microscopic disorder, so measurements of a single 
ring over a wide range of B can be interpreted as measurements over a large number of 



lithographically identical rings. This allows us to use the eight physically distinct rings to 
measure a much larger number of effectively independent rings. As described below, the large 
number of effective samples is crucial for making an accurate estimate of the higher cumulants 

of Pi?^ 

Based on these considerations, at large magnetic fields the persistent current is expected to 
take the form:^"* 

/(0„,0) = /W(0„)sin(2;za)/Oo) + /«(0„)cos(2;z^/Oo). [1] 

Theory makes three specific predictions concerning /"^^ and /"^ (the quadrature amplitudes of the 
AB oscillations). The first is that they are stochastic functions of Om characterized by a 
correlation function: 

(/W((D„)/<^>((D„ + AO^)) = {l'-\^J&\^^ + A^j) = (/2)C(A(D„ /OJ [2] 

that decays rapidly for AOm » Oc, where Oc is the correlation scale, which is typically a few 
times Oq. Both Oc and the normalized correlation function < C(x) < 1 have been calculated in 
Ref.[24]. 

The second prediction is that the distribution of these quadrature amplitudes = = Pj 

is Gaussian in the limit (g) ^ .^'^'^'^ For finite but large (g) it is predicted*' ** '^ that the 

normalized cumulant (defined below) of the persistent current k„ ~ g^'" (note that some k„ are 
suppressed still further by a large magnetic field^). In our samples g ~ 10"*, so the predicted 
deviations from Gaussianity are well below our present sensitivity. 

Lastly, correlations between /"^^ and /"^ are predicted to be absent. 

To test these three predictions, we first use the 1(B) data from one sample (ring #6, see SM 
for the full data sets and ring parameters) to determine the normalized autocorrelation of the PC, 

{l(B)I(B + AB)) / (^I^Y The result is plotted in Fig. 3, and shows AB oscillations whose envelope 

initially decays on a field scale ~ <E»c/27rrw, in qualitative agreement with the discussion above. 
After this initial decay the envelope does not approach zero, but instead undergoes apparently 



random fluctuations. These fluctuations are related to the flnite size of the data set, and are 
discussed further below. 

We can make a more quantitative comparison with theory by fitting the autocorrelation 
data in Fig. 3 to the prediction^"* that it should consist of AB oscillations whose envelope is given 
by C(AB/5c), where 5c = <^JlTtrw. The resulting fit is shown as the red line in Fig. 3. The fit 
parameters are S^''^ = 37 mT and = 25 mT (here the superscript denotes the ring #), in good 

agreement with the dimensions of the ring. The autocorrelation data from the other seven rings 
showed comparable agreement with the theoretical prediction, although the fitted values of Be 
varied from ring to ring (all values of the B^'^ are given in the SM). This analysis provides us 
with two useful results. The first result is the agreement between the measured and predicted 
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form of C(x), which justifies our use of the analytic expression for C(x) in the analyses below. 
The second result is the determination of the correlation field B^'^ for each ring, which will also 
be used below. 

To determine Pj from our measurements, we begin by applying the Hilbert transform to the 
1(B) data from each ring. This provides the quadrature amplitudes and f \B), as shown in 

Fig. 2(c). It is then straightforward to compute the cumulants of /"^^ and e.g., from their 
moments. Since there is no physical distinction between /"^^ and /"^ when Om » Oo, we consider 

the average of their cumulants: = ^ (((^'^^0 //'''((('^*^ // • account for variations 



between the rings (e.g., of D, r, and 7), the contribution to ^(^"// from each ring is normalized 

by the variance ^(^^)^ of that ring, giving the normalized cumulant = ((^ ))^((^^)) 
entire data set. 

The first several k„ are plotted as blue circles in Fig. 4(a). The prediction that P/is Gaussian 
(corresponding to /c„ = for all n > 3) is indicated by the black circles in Fig. 4(a). The data 
appear qualitatively consistent with a Gaussian distribution; however, to make a meaningful 
comparison between experiment and theory we estimate the uncertainty in these values. The two 
most important sources of uncertainty in the measurements of a:„ are the finite SNR of the 1(B) 
data and the finite size of the data set from which the Kn are calculated. We estimate the impact 
of the former by applying standard error-propagation techniques to the known uncertainty in 




each 1(B) measurement. This procedure is described in the SM, and leads to the blue error bars in 
Fig. 4(a). 

The finite size of the data set leads to a statistical uncertainty 6/c„ in the estimate of each k„. 
If the data sets for and ^'\B) each consisted of uncorrelated measurements, then values for 
the dxn could be found in standard statistics references. However it is clear from Fig. 2(c) and 
Fig. 3 that each quadrature of the AB oscillations contains correlations that are characterized by 
the function C(x). In this case the 6a;„ depend upon the form of C(x) and 5span, the range of B over 
which the PC is measured:^^ 
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c„ = j {c(x)rdx 



[3] 



(where the dependence of Be, B^pan, and c„ upon the ring # has been supressed). The black error 
bars in Fig. 4(a) correspond to the 6/c„ calculated from Eq. 3. Since 6/c„ < 1 only for n < 7, we 
plot the results only up to /ce. 

In addition to calculating the cumulants of Pi from our data, we can also plot the 
measured Pi in the form of a histogram of the /"^\fi) and I^'\B) data. To reduce oversampling 
artifacts in this histogram, we first bin the /"^\5) and f \B) data into a smaller data set. We 
choose the size of this smaller data set to correspond to the number of effectively independent 
data points in the entire data set^^ (i.e., from all eight physically distinct rings) 

8 

N^ff = 25^5^;^ / c^'>S<'' -412 (the factor of 2 in this expression arises from the two quadratures). 

(=1 

Here we have used n = 2 somewhat arbitrarily, but we note that A^eff depends only weakly upon 
the choice of n. As in the calculations of the the and /* \S) data from each different 

ring are normalized by their own variance to account for differences among the rings. The 
histogram of the resulting data set is shown in Fig. 4(b), along with the no-free-parameter 
prediction that this histogram should be Gaussian with zero mean and unit variance. 

We can apply these results to understand the magnitude of the apparently random 
fluctuations of the autocorrelation data in Fig. 3. At large AB, the data is expected to be 
uncorrelated (i.e., C(x) approaches zero for large x). However the standard error of the 



autocorrelation of a data set consisting of N independent samples is SC = l/ ^/]v . The data set for 
ring #6 contains =44 independent samples; thus at large AB the envelope of the 
autocorrelation should have a typical value « 0.15. This value is indicated by the dashed lines in 
Fig. 3, and is in agreement with the data. 

Lastly, we test our data for correlations between the quadrature amplitudes. From the 
f'^\B) and f\B) data it is straightforward to calculate the experimental value of 

cov(/'^\/'"^) 
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which is predicted to be zero. From our data we find p+. = 0.02 + 0.05, where the uncertainty 
arises from the finite size of the data set. This result is consistent with the prediction that 
and I^ \B) should be independent of each other. 

In conclusion, we have measured the distribution of mesoscopic fluctuations in metal 
samples. The form of this distribution is found to be consistent with theory, although more 
sensitive measurements would be required to detect the predicted non-Gaussian fluctuations. 
This result was obtained by measuring samples that each consist of a single coherent volume 
without any leads, and highlights the unique capability of persistent current to probe mesoscopic 
phenomena that cannot easily be accessed by traditional techniques. 

These results also test the independent-electron theory of persistent current (which has 
been a topic of controversy for more than two decades) in a novel regime: the persistent current's 
higher-order statistical properties. This opens the possibility of using persistent current to study 
other aspects of mesoscopic fluctuations such as multi-point correlation functions, which are 
predicted to show the direct precursors to Anderson localization even in the metallic regime.^ 

Lastly, the improved sensitivity of these experiments should enable the measurement of 
persistent current at much lower magnetic fields, where electron-electron interactions are 
expected to dominate the persistent current. ^^'^^ This is a regime in which many of the basic 
aspects of persistent current remain poorly understood. 
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Wohlmann, Leonid Glazman, Rob Ilic, Yoseph Imry, Felix von Oppen, and Douglas Stone for 
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Figure 1. Scanning electron micrograph of a cantilever with a single ring similar to the ones 
used in the experiment. Inset: magnified view of the ring. These images show the cantilever and 
ring prior to their release from the underlying SiOi layer. 
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Figure 2. Typical measurements of the persistent current, (a) A small section of 1(B), the 
persistent current as a function the applied magnetic field. The oscillations are due to the 
Aharonov-Bohm effect, (b) Red curve: 1(B) over a broader range of magnetic field. Blue curve: 
the envelope of 1(B). (c) The quadrature amplitudes l'^'^^ (black) and (pink). The envelope in 
(b) and the quadratures in (c) were extracted by applying the Hilbert transform to the red trace in 
(b), as described in the SM. 
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Figure 3. Autocorrelation of the persistent current. The blue curve is the normalized 
autocorrelation of the 7(5) data from ring #6, while the red curve is a fit to theory. Only data with 
AB < 0.3 was used for the fit. The expected error in the autocorrelation (due to the finite size of 
the data set) is indicated by the dashed horizontal lines. Similar behavior was observed in each of 
the seven other rings. 
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Figure 4. The measured distribution of persistent current, (a) Cumulants of the persistent current 
distribution, calculated from the combined data for all eight rings. Blue circles are the measured 
cumulants, and black circles are the cumulants expected for a Gaussian distribution. Black error 
bars: statistical uncertainty from the finite sample size. Blue error bars: statistical uncertainty 
from the finite signal-to-noise ratio, (b) Histogram of the observed persistent currents. To reduce 
oversampling in the histogram, the data is binned to give approximately independent data points. 
The solid hne is the no-free-parameters prediction. 
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1 Transport measurements 

The main paper compares our measurements of the persistent current in normal metal rings 
to predictions of a theoretical model based on diffusive, non-interacting electrons [3]. This 
theory depends explicitly upon two parameters of the rings: their diffusion coefficient D and 
their spin-orbit scattering length Lso- This theory also implictly assumes that the rings' 
phase coherence length {L^) is much greater than the rings' circumference L. The persistent 
current data presented in the main paper can be compared with this theory by taking D 
and Lso as fitting parameters (and assuming ^ L), but in order to constrain this 
comparison more tightly, we have used transport measurements to directly measure D, Lso, 
and L^. Transport measurements similar to those described in Ref. [2] were performed on an 
aluminum wire codeposited with the rings studied in this article. The wire was deposited on 
the same wafer as the rings. It had a length L = 255 /im, a width w = 105 nm (determined 
from SEM images) and a thickness t = 115 nm (determined from AFM measurements). 

Measurements of the wire's resistance were performed using a bridge circuit similar to 
that of Ref. [6]: a three terminal arrangement was used here as one of the four original 
leads was unintentionally blown out as the sample cooled to 4.2 K. One side of the sample 
was connected to the bridge and the voltage probe through separate leads, while the other 
side of the sample was connected to ground and the bridge through the same lead. The 
sample resistance was distinguished from the lead resistance by measuring the resistance of 
the sample as it transitioned from the superconducting to the normal state as the magnetic 
field was swept at a temperature below the wire's transition temperature Tc- A measurement 
of this transition is shown in Fig. [T^. 

The diffusion constant D was then calculated using the Einstein relation {p~^ = e^gD), 
with e being the electron charge and g the electron density of states per unit volume at the 
Fermi level for aluminum [7|. The resistance was measured to be 288 Jl. Based on the film 
dimensions, we infer a resistivity p = 1.36 x 10~® it 0.1 Om. This corresponds to a diffusion 
constant Dp = 0.02 it 0.0015 m^/s. We label Dp the diffusion constant extracted from this 
measurement of p. 

We also performed magnetoresistance measurements on the same wire at temperatures 
above Tc in order to extract and Lso- The magnetoresistance in the aluminum wire 
has contributions from weak localization as well as Maki- Thompson superconducting fiuc- 
tuations. The analytical form for the weak-localization correction to the resistance of a 1-D 

*To whom correspondence should be addressed; E-mail: manuel.castellanosbeltran@nist.gov. 
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wire in a magnetic field is given by 



where l/Lg = 1/-^^^ + 4/3-L|.q, is the electron phase coherence and Lso is the spin-orbit 
coupling length. The functions fi{B,Bl) and b{l) are defined in Ref. [2]. The Maki- 
Thompson contibrution is given by 

-n^]h{BML4.)) (2) 



R \Tc, 

where (3{t) is a function that diverges logarithmically when t — >■ 1 (see Ref. |2j and references 
therein). Equation |2] is only valid provided 

« ^ In (T/n) . (3) 

Fits were done using the sum of both equations [T] and |2| 

dR dR^^ dR"'^ 

— = \ (4) 

R R R 

Magnetoresistance measurements were made at a series of temperatures above Tc between 
1.8 and 12 K. Lso was obtained from fits of the magnetoresistance to Eq. |4] measured at 
the highest temperatures (T = 9 — 12 K), where Lgo was measured to be 1.54 it 0.06 ^mj^ 
For fits of the lower temperature data, Lgo was fixed at this value, so was the only 
fitting parameter. Fits to the magnetoresistace data using Eq. [4] are shown in Fig. [ija. The 
constraint set by Eq. [3] upon the validity of Eq. [2] sets the limit over which the fits were 
performed. 

The temperature dependence of L,^ can be explained based on the different contributions 
to the inelastic collisions of the electrons. Theory predicts that is limited by inelastic 
collisions with other electrons through the screened Coulomb interactions (ree), with phonons 
(rgp) and with extrinsic sources such as magnetic impurites. The latter should be negligible 
for the high purity aluminum source used and since no magnetic impurity has been observed 
to behave as a localized moment when dissolved in aluminum The fitted values of 
as a function of temperature are shown in Fig. The function used to fit to the data is 

= sJDt^, where 

= AepT^" + AeeT^I\ (5) 

A(,p and j4ee are fit parameters. The first term of Eq. |5] corresponds to the electron-phonon 
scattering rate, T~p = A^pT^- From our fit we find that Aep = 2.0 ± 0.2 x lO^s'^K-^ which 
is within a factor of two of previoulsly measured electron-phonon coefficients for comparable 
Aluminum wires jS] . 

The second term of Eq. |5] corresponds to electron-electron scattering. The wire (and 
rings) studied here had a width and thickness smaller than L^. Therefore, the quasi-lD 

'^The error bar in Lso was obtained from the following analysis: each of the fits to the four traces at the 
tempratures T = 9, 10, 11 and 12 K provided the same Lso within the fit error. The errorbars from the 
obtained values of Lso at each temperature are based on the goodness of the fits. Then we do a weighted 



average of the obtained values of Lso and the error is correspondingly a — l/y^S^^jT/a 
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Figure 1: (a) Resistance measurement of the aluminum wire as the magnetic field is swept 
through the critical field Be at a temperature T = 300 mK. The measured resistance in the 
normal state was 288 0,. Due to the loss of one the leads, in series with the sample, we 
measured a lead resistance of 22 0. This finite resistance is observed even when at \B\ < Be 
and T < Tc, and it has been subtracted from the trace, (b) Magnetoresistance data for T = 
2 K (o), 4.1 K (x) and 11 K (.) and fits (solid lines). The trace for 11 K has been multiplied 
by 20. At low temperatures, superconducting fluctuations are the dominant effect observed 
in the magnetoresistance measurements. The weak localization becomes the dominant effect 
at T > 8 K. (c) Phase coherence length L,^ extracted from fitting the magnetoresistance 
data as a function of temperature (points) and fit (solid line) to the expected coherence 
length assuming both electron-electron and electron-phonon contributions. Inset: Log plot 
showing Tg-p and Tge contribution as dashed-green and dashed-black lines respectively. 
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prediction for electron-electron interaction applies and the expression for Tee according to 
Ref. fa] is 



w 



2/3 



rp2/3 



(6) 



The fitted valued for 74ee is 0.65 ±0.03 X lO^s The expected theoretical value based 

on Eq. [ojand a diffusion constant of L> = 0.020 mVs is Aee = 0.15 x lO^s-^R-^/s. Although 
there is a large discrepancy betwen the expected and fitted value of Aee, the primary purpose 
of these transport measurements is to show that is greater than the circumference of the 
rings L = 1.8 and 2.8 //m; therefore we did not look further into this disagreement]^ 

The value for the diffusion constant extracted from the resistivity measurement can then 
be compared with the diffusion constant extracted from the persistent current measurement 
(Fig. [2]). Using a similar analysis to that explained in Ref. [3] we extract a value of 
Dpc = 0.021 lb 0.002 m^/s. We label this value as Dpc to distinguish from Dp. Dpc and 
Dp differ by about 5%. This seems reasonable given differences between wire cross section 
and ring cross section plus statistical uncertainty in Dpc- 



2 Data Analysis 

We inferred the persistent current in each ring by measuring the resonance frequency of the 
cantilever into which the ring was integrated by driving the cantilever in a phase-locked loop 
[10] . The method used to infer the current is explained in the supplementary information 
of Ref. [2]. Here we briefly review this analysis. 

As it was shown in Ref. fH], the presence of magnetic field inside the metal of the 
ring leads to an aperiodic modulation of the persistent current oscillations and a change in 
the flux dependence of of the persistent current from the simpler case where there is only 
a pure Aharonov-Bohm flux (j) threading the ring. This change is the modification of the 
time-reversal relation from /((/)) = —I(—(p) to I(Bm,(p) = —I{—Bm, — </>) . As a result, the 
current is no longer odd in the Aharonov-Bohm flux cf) and it takes the more general form 
at fixed Bm 

m = E It/,e «i-(2-f ^) + lH„e ^°«(2-f ^) (7) 
P 

where the variables iti and 17, are stochastic variables that vary with the magnetic field 

n/pe n/pe ° 

B^ as well as with microscopic disorder, and $0 = ^/e is a flux quantum. Determining the 
distribution of ^^i^^ and l^j^^ is the central point of the main paper. 

We monitor the the cantilever frequency as we vary the magnetic field. The cantilever's 
deflection leads to a rotation of the sample, responsible for the coupling between the per- 
sistent current and the cantilever. When 7^ 7r/2 (where Q is the angle between the plane 
of the ring and the applied magnetic field B\ the frequency change is dominated by the 
following term: 

= -4^cos(^)Vjif2.rp^^^|^^V (8) 

, , BA^avlQ. _ , BAsmQ. 

hipe cos(2^p^— ) - 4/^^ sin(2^p^— . 



^ In addition, it has been pointed out that close to the superconducting transition, the electron-electron 
inelastic scattering rate can be modified due to superconducting fluctuations [8 . These fluctuations may 
alter the numerical coefficient Aee- 
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Figure 2: (a) Frequency shift versus magnetic field for a ring of radius r = 296 nm at 
different temperatures. The amphtude of frequency oscillations decreases with increasing 
temperature as expected. The amplitude of the frequency oscillations is converterted into 
current oscillations amplitude and plotted as a function of temperature in (b). Data (circles) 
and a fit (solid line) based on analysis similar to that explained in Ref. [3] are plotted. From 
the fit we estimate a diffusion constant of D = 0.021 it 0.02. 
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where A is the area of the ring, /q is the natural resonance frequency of the cantilever, 7^ 
is the ratio between the slope of the cantilever and the factor xup/L for the flexural mode 
m. For m = {1, 2}, 7^ = {1.377, 4.788}. k is the cantilever spring constant, L is the length 
of the cantilever, xtip is the amplitude of oscillation at the tip of the cantilever, and Ji(x) 
is the first Bessel function of the first kind. 

At the temperatures of our experiment the current is dominated by the first harmonic, 
p = 1. Thus, for the analysis of the data shown in the main paper we use what in the 
supplementary online information (SOI) of Ref. |2] is refered to as method B: we assume 
that the signal only has the p = 1 component in Eq. [8] and that the argument of the Bessel 
function varies only weakly over a given data set. In that case, the change in frequency is 
essentially the derivative of the persistent current (again ignoring all terms for p > 1). 



dl , 27r^sin6' 
of 



dB ' 



ABv^yn T fn AB COs{9) JrnXtip 

cos(t^) Ji I zirp 



kLxtip \ ^0 L 



-1 



In order to estimate the current this quantity can be numerically integrated. However, since 
we are interested in the statisitics of the variables Ij^^^, we perform the analysis on /'defined 
as 

/'(B) = (9) 

zttA sni & oB 

The normalization is such that the oscillations of I'{B) have the same amplitude as those 
in I{B). The reason we use /' is that the numerical integration can introduce unwanted 
correlations in the values of the variables If^, . 

h/e 

2.1 Drift removal 

We measured cantilevers similar to the ones shown in Fig. |3] We fabricated long (/ ~ 400/im) 
and short (I ~ 100 /xm) cantilevers, with three different widths (w = 20, 40 and 60 ^um). The 
signal 6f was larger for the shortest cantilevers. However, their noise peformance was very 
poor, presenting frequency noise with a power dependence ~ 1//" with n ~ 2, preventing 
our measurements from achieving the thermal noise limit. However, for the long cantilevers 
the noise performance of the frequency measurement did reach the thermal noise limit (see 



Sect. 2.2.1). 
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Figure 3: Optical micrograph of cantilevers similar to those used for the measurement of 
persistent current in the main text. The ring sample on each cantilever is not visible. 



Raw data of the frequency measurements of two kind of cantilevers are shown in figures 
[4] and [5] In both figures we can see a drift in the cantilever's resonance frequency with time, 
but with very different characterisitics. In Fig. |4] the Aharonov-Bohm (AB) oscillations 
can be easily distinguished on top of a mostly magnetic- field dependent drift. This drift is 
removed using MATLAB's local regression algorithm for the LOWESS (Locally Weighted 
Scatterplot Smoothing) routine with a first degree polynomial. We have observed that if we 
choose the window of the LOWESS routine to be the equivalent of 5 AB oscillations, the 
peak of our signal in its Fourier transform is unchanged by the drift removal process. 

For a thermally limited frequency measurement, shorter cantilevers should have a greater 
signal-to-noise ratio, but we found that, in practice, their frequency noise was significantly 
worse than the longer cantilevers. Thus, for shorter cantilevers, the frequency measurements 
were not thermally limited. This is obvious in figures and b where it is not possible to 
discern the AB oscillations at first sight. In Fig. ^ and d we show the averaged data and 
their spectral densities with and without LOWESS drift removal in order to show that our 
procedure does not convert a noise shoulder into a peak. 
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Figure 4: (a) Cantilever frequency (blue left axis) as a function of time as the magnetic field 
(green right axis) is varied for sample 6. On top of a magnetic-field-dependent background, 
the oscillations due to the persistent current are obvious and the drift in the cantilever 
frequency is easily removed by using local regression. After substracting the drift calculated 
using LOWESS (red trace) we obtain the oscillations shown in (b). 
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Figure 5: (a) Cantilever frequency (blue left axis) as the magnetic field (green right axis) 
is varied for sample 2. In this case, the time dependent background is so large, it is hard 
to distinguish the AB oscillations even after removing the drift as shown in (b). (c) The 
AB oscillations can be easily distinguished after averaging all the data with the drift (red) 
and with the drift substracted (blue). This proves that the observed oscillations are not 
an artifact of the procedure used to remove the drift, (d) Spectral density of (c) showing 
that the procedure used to remove the drift does not affect the amplitude of the observed 
oscillations. 



2.2 Hilbert Transform 

The Hilbert transform of a function u{t) is defined as: 

p.v. — (g) u(t) = H[u(t)] = u(t) 
where p.v. is the Cauchy principal value and indicates convolution. The Bedrosian 
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theorem states that the Hilbert transform acting on the product of two functions x{t) = 
A{t) X h{t) can be written as 

H[A{t)h{t)] = A{t)H[h{t)] (10) 

if the Fourier spectra for A{t) and h{t) are disjoint in frequency space and if the spectrum 
for h{t) is concentrated at higher frequencies than those of A(t). For example, if h(t) is a 
periodic function h{t) = cos (j){t), h{t) = sm(j){t). In this regime, then the following relations 
hold dl]: 

x{t) = A{t) cos (p{t) (11) 



A{t) = ahs{x{t) + ix{t)) (12) 

cp{t) = arg(x(t) + ix{t)) (13) 

For the work described in this paper, the Hilbert transform is used to analyze two 
separate data sets. In the first case, the Hilbert transform is applied to raw interferometer 
data (i.e. cantilever position vs. time). This technique allows us to generate densely 
sampled frequency vs. time traces, which are used to diagnose the sources of noise in the 



measurement. We discuss this technique in section 2.2. 1| In the second case, the Hilbert 



transform is used to extract the persistent current quadrature amplitudes I^^^ from the 



persistent current data. This is discussed in section 2.2.2 



2.2.1 Phase Noise Analysis 

As explained in Ref. ||2], the persistent current measurement is at its core a frequency (/) 
measurement of a driven cantilever. In order to better understand the specific noise sources 
of the measurement system, it is beneficial to have access to the noise spectral density Sf{uj) 
of raw frequency data f{t) 

Sfiu:) = 2 lim 
where Ft{uj) is the windowed fourier transform of f{t) 

T/2 

iujt' 



Ft{u:)= I dt'f{t)e 



-T/2 

The factor of 2 comes from the fact that we are only considering single-sided spectral densities 
(thus we only consider a; > values). For the work described in the main text, we utitlize 
a technique involving the Hilbert transform to measure the frequency f{t). The Hilbert 
transform is applied to cantilever position data obtained by an optical-fiber interferometer. 

If the interferometer data contains more than one frequency component, we can define 
an "instantaneous frequency" f{t) = fo + fnoise{t)- Here fnoise{t) is a stochastic variable with 
zero mean and /o is a constant. We will consider the limit in which the Fourier transform 
of the cantilever position is sharply peaked at /o and also /o ^ A/, where /o is the center 



frequency and A/ is the width of the peak. In this limit, the Bedrosian theorem (Eq. 10) 



holds, and we can use the Hilbert transform to calculate the cantilever phase vs. time with 



Eq. 13 
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We convert the phase versus time (/>(t) of the cantilever's motion to an instantaneous 
frequency f{t) via the relationship: 



2,m = i (") 

The frequency noise spectrum is related to the phase noise spectrum by Sj = {(jj/2'k)'^S^ by 



Eq. (14). We can then compare the measured Sj to theoretical predictions for the frequency 
noise of a driven limit-cycle oscillator subject to a white force noise and a white displacement 
noise in the detection 1101 



thermal (w) 



tip 



25rV^\^ ^^^^ 

•J/, detector (Wj = 2 \o~ I ^ > 

^tip ^^'^^ 

where k is the resonator spring constant, Q is the resonator quality factor, Xtip is the 
resonator displacement amplitude, and S'™'' is the displacement noise of the detector. The 
imprecision in the measurement of the cantilever's motion arises from fluctuations both 
in the laser used to monitor the cantilever and in the detector used to measure the laser 
signal. The dominant noise source in our case is the electronic noise from the photoreciver. 
In order to convert the voltage at the output of the photodetector into cantilever motion, we 
measure this voltage when the cantilever is only excited by a thermal force. The magnitude 
of the cantilever's thermal motion can be computed from the equipartition theorem which 
states that at thermal equilibrium 

k{x^) = kBT (16) 

The spectrum of the voltage Sy at the output of the photodectector has the shape of 
lorentzian curve on top of an offset (inset of Figure |6]). This offset consists of the mea- 
surement imprecision S™*^. By measuring Sy at different temperatures (Figure [6|, we can 
properly convert the volts at the output of the detector into displacement of the cantilever 
and calculate the imprecision ST" of our detector. 

In these measurements, for some cantilevers, we notice a strong deviation from the 
prediction at low frequencies. This low-frequency behavior seems to correlate with the 
amount of mechanical nonlinearity in the cantilever (determined, e.g., by a non-Lorentzian 
resonance for large Xtip). For samples 3-8, there was very little nonlinearity and the frequency 
noise behaves as expected (Fig. [7]). For samples 1 and 2, we noticed a large amount of 
mechanical nonlinearity and also a large amount of low- frequency noise (Fig. [s]). 

2.2.2 Calculating the Persistent Current Amplitude 

The goal of the work in the main text is to measure the statistics of the quadrature am- 
plitudes of the AB oscillations I^j^- This amplitude slowly changes due to the magnetic 
field that penetrates the metal of the ring 0. Since our data is dominated by the first 
harmonic of Eq. [7] (no higher harmonics are visible in the data) the normalized derivative 
of the current /' (Eq. [9]) takes the form: 

= /+,cos(2vr^,^)-4-/^sin(2^p^) (17) 
= h/e{B)cosi27rp^— + a{B)\ (18) 
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Figure 6: Thermal motion of the cantilever in units of temperature (Eq. 16 ) as a function of 
the fridge temperature. This motion is measured as a voltage at the output of the photode- 
tector and converted into displacement units (inset). The spectral density of the voltage has 
two constributions: a white background that consists of the measurement imprecision Sx ^ 
coming from the voltage noise of our photodetector (dashed blue line in the inset) , and a 
lorentzian curve that consists of the thermal motion of the cantilever. Measuring the detec- 
tor's voltage fluctuations at different temperatures allows us to calibrate the intereferometer 
response. 
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Figure 7: Frequency noise for sample 4. The solid lines represent the frequency noise 
spectrum data obtained using the Hilbert transform technique described above at various 
oscillator drive strengths. The dashed lines are the theoretical noise curves calculated with 
Eq. 15 (no fit parameters). There are two major contributions to the noise: thermal 
fluctuations driving the cantilever, and the intrinsic detector noise. The detector noise was 
extracted from the actual interferometer noise. For this cantilever, Q = 5700, T = 320 mK, 
/o = 966 Hz and ki = 4.0 x 10"^ N/m. 
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Figure 8: Frequency noise for sample 2. The solid lines represent the frequency noise 
spectrum data obtained using the Hilbert transform technique described above at various 
oscillator drive strengths. The dashed lines are the theoretical noise curves calculated with 
Eq. 15 (no fit parameters). This sample exhibited strongly nonlinear behavior and its 
low-frequency noise deviates from the theoretical curves. For this cantilever, Q = 35000, 
T = 320 mK, /o = 11827 Hz and ki = 1.7 x 10"^ N/m. 
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Figure 9: (a) Hilbert transform for the data from sample 6. From this trace we obtain 
Ih/e{B) and a (not shown), (b) Using the data from (a) we can infer the values for /^^^ and 



where Ifi/^{B) is the amplitude of the current and cos \^2Trp — a{B)j represents the 

AB oscillations. 

The correlation field Be sets the field scale over which Ii^jf,{B) changes. In the limit 
that Be is "large enough" [Be 3> ^q/A)^ Ii^ie{B) is a slowly-varying function compared 



to cos(27rp^^^|^^ - a) and the Bedr osian theorem pO^ holds. The AB frequency is also 
sharply peaked in Fourier space, with a well-defined irequency given by the dimensions of 



the ring. Thus, comparing (12) and (17), we can obtain If^/f,{B) and a using the Hilbert 
transform. This also allows us to determine the quadrature variables 1^^^ = Ih/^ cos(a) and 
^h/e ~ -^/i/e sin(a). This technique is illustrated in Fig. joj 

2.3 Definition of the cumulants and statistical moments 

For the purpose of this manuscript, it is convenient to have a consistent definition of the var- 
ious moments and cumulants. The raw moments //^ and central moments fir of a stochastic 
variable x are defined as 

f^r = {^1 (19) 
Hr = {{x-n'iY) 

The cumulants ((/")) are most easily defined in terms of the central moments: 

((^^» = 



15 



{{I")) = (20) 
= ^2 - 15^4^2 - 10/ii + 30//! 

2.4 Finite sample statistics 

In Ref. [3j it was mentioned that in the presence of an additional large in-plane magnetic 
field Bm penetrating the metal ring, the effective disorder of the ring changes, implying that 
averaging over magnetic field is equivalent to an ensemble average [3]. However, our finite 
magnetic field range means that in practice, we have a finite number of realizations from this 
ensemble. In order to estimate the statistical uncertainty in our estimates of the cumulants 
due to this finite sample size, we use the results of Ref. [5J. 

The statistical uncertainty of the cumulants (((/"))) due to a finite sample size can be 
expressed in terms of the normalized correlation function C{x) defined as 

^,'B-B'\ _ {mi{B')) 



Be J {P) 

where C decays from C(0) = 1 to C{oo) = 0, and Be is the correlation field, which sets a 
rough order of magnitude over which the persistent current is correlated. Expressions for C 
are provided in Ref. The typical value of a cumulant calculated from a finite ensemble 
is given by 

where = n\ /^[C(a;)]"'da;, and -Bgpan is the total magnetic field range over which data is 
taken. The theoretical prediction for the higher order cumulants of the persistent current 
quadratures is 

gn Z 

for n > 2, where the typical current is defined as /typ = \/ (/^) [1]. In our case g ~ 10^ and 
thus we expect a gaussian distribution for I^j^ and the statistical error is thus given by 

H{n) = Janj^{pr (21) 

For an estimate of this term, we can consider the case for a spinless electron at T = 
and large Bf^ (so that we can ignore the Cooperon contribution) and only include the first 
harmonic of the PC. In this case 



C(x) = (^1 + |x| + yj e-l"l (22) 

For this simplified case, the first 4 coefficients are 02 = 7, 03 = 16.61, 04 = 56.48, and 
05 = 249.55. Experimental values of the cumulants of the persistent currents are considered 



sound only if they comfortably exceed the systematic error ycin ^ " (/^)" (Eq. 21). The 
actual correlation function for our case is more complicated than Eq. [22] and is a function 
of T and Lso- It can be numerically implemented using Ref. [3J. A fit to the correlation 
function C{B / Bc,T, Lso) of our data is shown in Fig. 10^ 



^The correlation function is a also a function of the Zeeman energy, but at the large magnetic fields used 
in this experiment, this dependence is negligible. 
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Figure 10: Persistent current autocorrelation for sample 6. Blue represents the autocorrela- 



tion of the data shown in Fig. 20 and red is the fit using Ref. |3]. The estimated correlation 
field in this case is Be = 32 mT. In order to estimate the standard error in the autocorrelation 
(dashed horizontal lines), we assume that this is approximatly given by the standard error 
for the correlation coefficient l/^/N^ ~ 0.15. As it can be seen in the figure, fluctuations 
of this order are present at large AB ^ Be, where the autocorrelation of an infinitely large 
data set would be expected to vanish. 



We can also define an effective number of samples based on the expected statistical 
uncertainty for the cumulants of an ensemble of A'eff samples |12) : 



6 an) 



/n!((l2))" 



By comparing the latter with Eq. 21 , we can then define A^eff as 

n\B, 



-'span 



an Be 



(23) 



2.5 Finite signal to noise 

Another source of error in our estimate of the cumulants is the finite sensitivity of our 
setup. In order to estimate this error, we have to calcualte the appropriate expression for 
error propagation in our analysis. Our case is somewhat confusing, since Xi itself is an 
stochastic variable, so really we have Xi^meas = Xi + Sxi. First let's consider the general 
equation for the error propagation for a general function f{xi, X2-.-,X]\f). If Xi is known only 
to within an error 6xi then the error in / is given by 

i ^ ' 

Let's first consider the case of the central moment as cumulants are easily definined in terms 
of them: ^ 

f{Xi) = = — ^ (xi - /i'l)'' 
i 
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Then 



N 



f^2r-2 {6^X) (25) 



There are two assumptions used for Eq. [25] First, we assume that the noise 5xi from the 
different magnetic field points (labeled with the index i) are uncorrelated. This does not 
mean that the different uncorrelated. The second assumption is that 6 Xi can be 

replaced by an average Although not completely accurate, since the sensitivity of our 
measurement is lower at lower magnetic fields, we did compensate for the loss of sensitivity 
by averaging longer at low magnetic fields. 

The expressions for the second and third cumulants are the same as for the second and 
third central moments. Thus, Eq. [25] gives their measurement uncertainty. For the forth, 
fifth and sixth cumulant we can use the following heuristic approach |12) . We start with the 



definition of the cumulants (Eq. 20), from which the following expressions are derived: 



6 {{l^}) = - 15/146^12 - 15^2(5^4 - 20//35/i3 + 90^16^2 
Then, for example, for the fourth cumulant: 

(52 = ,5^4 + 36fil6^H2 - 12/X2 {5^14S^12) (26) 

The covariance {5 1x^^8112) is not necessarily zero and it can be derived using a similar approach 



as the one used to derive Eq. 25 



Similar expressions to Eq. |26|for ((J^)) and ((/*')) can be derived. 



2.6 Data 

Following, we show a table with all the estimated cumulants and estimated errors for the cu- 
muluants of all 8 samples. The cumulants of both /"*" and I~ have been combined for each of 
the rings. In order to account for variations between rings the different cumulants are normal- 
ized by the variance (^(^I^'^y, thus we define a normalized cumulant k„ = ((/")) / ((l^))"''^. 
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^QTYIT^ £i -H- 
OOjlliL/iC -j-f- 


/ 1 lllli J 






//Cr±\2\ \ / A2^ 
\\ // 


^3 


1 


296 


16 ±1 


63 


0.46 X (1± 0.18 ±0.02) 


0.065 ±0.28 ±0.05 


2 


296 


25 ± 1.5 


43 


0.37 X (1± 0.22 ±0.02) 


0.23 ±0.34 ±0.04 


3 


448 


9± 1 


55 


0.0031 X (1 ± 0.2 ±0.03) 


0.21 ±0.31 ±0.08 


4 


448 


17 ± 1 


34 


0.0063 X (1 ± 0.24 ±0.02) 


-0.25 ±0.4 ±0.07 


5 


296 


22 ± 1 


56 


0.597 X (1 ± 0.19 ± 0.01) 


0.18 ± 0.29 ± 0.03 


6 


296 


32 ±2 


44 


0.67 X (1± 0.21 ±0.01) 


0.001 ±0.34 ±0.02 


7 


296 


20 ±1 


81 


0.26 X (1± 0.16 ±0.01) 


0.083 ±0.24 ±0.03 


8 


418 


11 ± 1 


36 


0.022 X (1 ±0.23 ±0.02) 


-0.05 ± .37 ±0.05 



Sample 7^ 


k;4 


K5 


Kg 


1 


-0.052 ±0.51 ±0.1 


-0.12 ±1.1 ±0.25 


-0.91 ±2.6 ±0.47 


2 


-0.0271 ±0.62 ±0.07 


-0.83 ±1.3 ±0.2 


-0.93 ±3.1 ±0.41 


3 


0.29 ±0.58 ±0.18 


-0.05 ± 1.23 ±0.5 


-0.84 ±2.9 ±1.13 


4 


0.56 ±0.75 ±0.16 


0.27 ±1.58 ±0.5 


-0.043 ± 3.7 ± 1.32 


5 


-0.25 ±0.54 ±0.04 


-0.17 ± 1.15 ±0.14 


0.36 ±2.7 ±0.28 


6 


-0.61 ±0.637 ±0.03 


-0.05 ± 1.35 ±0.09 


2.23 ±3.15 ±0.18 


7 


0.065 ± 0.45 ± 0.07 


0.51 ±0.95 ±0.23 


2.85 ±2.23 ±0.45 


8 


0.2 ±0.68 ±0.086 


0.12 ±1.4 ±0.23 


-2.1 ±3.35 ±0.58 



Table 1: Higher order moments of the quadrature amphtudes of the persistent current 
for all the measured samples. The first error indicates the standard error due to finite-size 
sample and the other represents uncertainty from noise in the frequency measurement. A'^2,eff 



indicates the effective number of samples for a given ring using Eq. 23 multiplied by 2 as 
we use both quadratures to calculate the cumulants. 



3 Measurement diagnostics 

We performed a set of diagnostic measurements similar to those described in Ref . j2] . Specif- 
ically, we measured the effects of the readout laser power on the measured frequency, we 
compared the extracted current using the two first modes of the cantilever motion and we 
checked the effects of cantilever oscillation amplitude. 



Figures 11 and 12 show the effect of varying the laser power upon the persistent current 
data measured on two different cantilevers for different incident laser powers. For different 
cantilevers, we observed two qualitatively different types of dependence upon the laser power. 
We believe this difference is due to the different widths of the cantilevers. The data of samples 
1-4 where taken with an incident power of 10 nW. However, for samples 5-7, the data was 
taken with an laser incident power of 3 nW as they presented a stronger power dependence. 
Sample 8 was taken with a laser power of 5 nW in a previous cooldown (Ref. [2]). 
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Sample 7^ 


r (nm) 


w (nm) 


t (nm) 


^cantilever (Z^^^) 


^cantilever (l^m) 


1 


296 


90 


115 


114 


40 


2 


296 


90 


115 


126 


40 


3 


448 


90 


115 


395 


40 


4 


448 


90 


115 


398 


40 


5 


296 


90 


115 


370 


20 


6 


296 


85 


115 


359 


20 


7 


296 


95 


115 


352 


20 


8 


418 


85 


90 


438 


60 



Table 2: Sample parameters. All the measurements were done at a 45-degree angle between 
the cantilever and ring plane and the magnetic field. All the cantilevers had a thickness of 
118 nm except for sample 8 which had a thickness of 340 nm. Sample 8 was measured in a 
separate cooldown. 



In Fig. 13 we demonstrate that the inferred current is the same whether the cantilever's 
first or second flexural mode is used, indicating that persistent current is independent of 



excitation frequency of the cantilever. Finally, in Fig. 14 we show that the cantilever's 
frequency depends upon the amplitude of its motion as would be expected if the persistent 
current remains in its equilibrium state (Eq. [8|. To generate the data shown in Fig. [l4j the 
resonant frequency of cantilever 5 was measured at two different magnetic fields (indicated 
in the inset figure with two arrows) as a function of cantilever amplitude Xtip. Then, the two 
measured 5f were substr acted in order to remove any kind of amplitude-dependent change 
in the cantilevers resonance frequency. The cantilever was excited in its first flexural mode. 
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Figure 11: Left panel: Change in frequency versus magnetic field for a series of laser powers 
incident on the cantilever. This data was taken with sample 6. The drive was the same for 
all the traces. Right panel: The mean amplitude of each of the traces from the left panel, 
plotted versus laser power. For the data of samples 5, 6 and 7 shown in the main paper, 
3 nW of laser power was used. 
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Figure 13: The derivative of the persistent current I' (derived from Eqs. [s] and [o]) ver- 
sus magnetic field measured when oscillating the cantilever at 1.232 kHz (the cantilever's 
first flexural resonance) and at 7.77 kHz (the cantilever's second flexural resonance). The 
persistent current does not appear to depend on the cantilever oscillation frequency. 
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Figure 12: Left panel: Change in frequency versus magnetic field for a series of laser powers 
incident on the cantilever. This data was taken with sample 2. The drive was the same for 
all the traces. Right panel: The mean amplitude of each of the traces from the left panel, 
plotted versus laser power. For the data of samples 1, 2, 3 and 4 shown in the main paper, 
10 nW of laser power was used. 
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Figure 14: Cantilever drive test on sample 5 where we measure the accuracy of Eq. [8] The 
cantilever amplitude is plotted on the x-axis in terms of the amplitude of the flux modulation 
4'ac/4'0 through the ring produced by the cantilever motion. Data points reperesent the 
difference in cantilever frequency shift for the two field values indicated in the inset. The 
solid curve is a fit using Eq. [s] (with p = 1) and r = 301 it 2 nm, consistent with the measured 
radius and linewidth of our ring (see table [2]). Inset: the arrows indicate two field values 
at which measurements of the cantilever frequency shift were performed as a function of 
cantilever amplitude. 
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4 Data of all the samples 



4.1 Magnetic field sweeps 



In this last section, we present figures 15 22 where we show the complete I' versus magnetic 
field traces which were analyzed in the main text. These traces were calculated, using method 
B of Ref. [2J from measurements of the cantilever frequency performed at the refrigerator's 
base temperature of 320 mK (365 mK for sample 8). 






B (Tesla) 



Figure 15: The derivative of the persistent current I' (derived from Eqs. [s] and [9]) versus 
magnetic field for sample 1 at T = 320 mK. The full sweep is separated into three contiguous 
panels for clarity. 
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B (Tesla) 

Figure 18: The derivative of the persistent current I' (derived from Eqs. [s] and [9| versus 
magnetic field for sample 4 at T = 320 mK. 
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Figure 19: The derivative of the persistent current /' (derived from Eqs. [s] and [9]) versus 
magnetic field for sample 5 at T = 320 mK. 
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Figure 20: The derivative of the persistent current /' (derived from Eqs. [s] and [9]) versus 
magnetic field for sample 6 at T = 320 mK. 
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